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Abstract. We present and analyze in detail an unknown theory of fcrroinag- 
netism developed by Ettore Majorana as early as the beginnings of 1930s, sub- 
stantially different in the methods employed from the well-known Heisenberg 
theory of 1928 (and from later formulations by Bloch and others). Similarly 
to this, however, it describes successfully the main features of ferromagnetism, 
although the key equation for the spontaneous mean magnetization and the 
expression for the Curie temperature are different from those deduced in the 
Heisenberg theory (and in the original phenomenological Weiss theory). The 
theory presented here contains also a peculiar prediction for the number of 
nearest neighbors required to realize ferromagnetism, which avoids the corre- 
sponding arbitrary assumption made by Heisenberg on the basis of known (at 
that time) experimental observations. Some applications of the theory (linear 
chain, triangular chain, etc.) are, as well, considered. 



1. Introduction 

The microscopic interpretation of the magnetic properties of bodies, related to the 
motion of electrons in atoms, appeared as early as the beginning of the XX century, 
when P. Langevin [T] provided a successful theory of diamagnetism, by evaluating 
the induced magnetization in presence of an applied magnetic field in the frame- 
work of the Lorentz theory of the electron. Subsequently the paramagnetism of 
materials was rightly related to the fact that atoms or molecules can have perma- 
nent dipole moments arising from the current loops generated by orbiting electrons 
(proportional to the total angular momentum), and thus acting as tiny magnetic 
shells. The correct interpretation of the ferromagnetic phenomena, instead, did not 
come out soon, and their understanding on an atomic basis was achieved only after 
the advent of quantum mechanics. 

The striking property to be explained was that the ferromagnetic substances 
(like iron, cobalt, nickel, etc.) develop a spontaneous magnetization only when 
cooled below a certain temperature T^, called the Curie point. Differently from 
an ordinary paramagnet, above this critical value the temperature variation of 
the magnetic susceptibility % of a ferromagnet was found empirically to follow the 
Curie- Weiss law, x = C / {T —Tc) (where C is a constant), whereas for a paramagnet 
the Curie law x — C jT holds. 

Early works by J. A. Ewing, K. Honda and J. Okubo, and otherfQ showed that 
many of the peculiarities exhibited by ferromagnetic bodies could be described by 
assuming a large potential energy between adjacent microscopic (atomic or molec- 
ular) magnets. However, a first important step towards the understanding of the 
phenomenon was performed by P. Weiss in 1907 [3], when he postulated that the 



^For a description of these and other earlier models see the review in Ref. [2]. 
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Figure 1. Graphical solution of the implicit equation llj, ob- 
tained from the intersection of the two curves M = Mq tanh x and 
M = (kT/ ^a)x for zero external magnetic field and different values 
of the temperature. A non-vanishing spontaneous magnetization 
(corresponding to the point P) exists only for T < Tc. 



atomic magnets tend to be brought into parallel orientation not only by an ap- 
plied, external magnetic field but also by an inner field Hi, called by Weiss the 
"molecular field", which is proportional to the magnetization M of the material: 
Hi = aM. By introducing such an idea of an inner field into the Langevin theory 
of the magnetism of classical dipoles, Weiss was in fact able to account for most 
of the phenomenology of ferromagnets [4J. If a dipole is brought from outside and 
placed in the molecular field, it will be aligned in the direction of this field (i.e. in 
the direction of the magnetization itself), giving rise to a spontaneous ordering of 
the elementary magnets. Such an ordering effect competes with random thermal 
effects, with energy of the order of kT, trying to flip the dipole away from the 
ordered state, but Weiss showed that (in the simple case that dipoles can align 
only in two opposite direction), even in the absence of an external magnetic field 
[He = 0), a non- vanishing magnetization may arise, satisfying the equation: 



(/i is the magnetic moment of the N atoms per unit volume of the body, and 
Mq — A^/i)0 In fact, apart from the trivial solution M = 0, the implicit equation ^ 
admits also another non-zero solution for the magnetization M , provided that T < 
Tc = Nfi^a/k, as it was seen graphically (see Fig. [T]). The magnetic susceptibility 
was as well calculated, obtaining the experimentally observed Curie- Weiss law (with 
C = N^'^/k), thus finally showing the good agreement of the phenomenological 
Weiss theory with the known properties of ferromagnetism. 

However, despite such a success, the question about the origin of the local field 
remained unanswered, and Weiss himself realized that it could not arise from the 
magnetic interactions of the magnetic moments of the atoms or molecules, these 



^Actually, in the original formulation, the Langevin function cosh 2: — 1/x appeared in Eq. JTJ 
instead of the function tanhx, which was considered only later by W. Heisenberg (see below). 
However, as it is easily recognizable (as done by Heisenberg himself), the physical interpretation 
of the phenomena does not change and no alteration in the reasoning is required, so that we prefer 
to use the mathematical form of Eq. |(TJ for later convenience. 



M = Ma tanh 



fj.aM 



(1) 



kT 
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being extremely small. In fact, the value of the proportionality coefficient a should 
be assumed as high as of the order of 10^, against the value of 47r/3 calculated under 
ordinary electromagnetic assumption [5]. On the other hand, classical electrostatic 
forces may lead to interactions of the right order of magnitude, but do not give a 
linear proportionality between the Weiss molecular field Hi and the magnetization 



The solution of this puzzle had to wait until the end of 1920s, when W. Heisen- 
berg [7] showed that the forces leading to ferromagnetism are purely quantum 
mechanical in nature, being due to electron exchange. The key idea was to show 
how strong non-magnetic forces between electrons that favor spin alignment (and, 
thus, the alignment of the elementary magnets) arise from the quantum mechani- 
cal exchange interactions, in the same way as they produce level splittings between 
electrons in singlet (para-) and triplet (ortho-) states of two-electron atoms: "the 
empirical phenomenon that ferromagentism presents is very similar to the situation 
we met earlier in the case of the helium atom." 0. Heisenberg, indeed, began in 
1926 [8J to explore how the so-called "exchange term" appears in a system of atoms 
with many electrons satisfying the Pauli exclusion principle. Exchange forces are 
electrostatic in nature but, because of the constraints imposed by the Pauli princi- 
ple, they are formally equivalent to a tremendously large coupling between spins, 
as seemingly required by the Weiss phenomenological theory. With a suitable as- 
sumption (a positive exchange integral), Heisenberg showed that the state in which 
spins are aligned in the same direction is energetically favorable, helping in the 
ordering of the spins in a certain direction, thus proving his model to be well fitted 
to describe ferromagnetism. 

This result was, actually, achieved by following the statistical model built up by 
E. Ising few years before [9]. In the Ising model (or, rather, the Lenz-Ising model 
[10]) of a ferromagnetic body, N elementary magnets of moment are assumed to be 
arranged into a regular lattice, each of them having only two possible orientations, 
as in the Weiss theory. The magnets experience only short-range interactions on 
each other, the interaction energy (assumed to be a constant) being non-vanishing 
and negative (positive) only for each pair of neighboring magnets of same (opposite) 
direction. This interaction tends, then, to make neighboring spins the same (an 
additional energy ±fiHe is present for an applied magnetic field He), so that in 
principle this model could be suitable for describing ferromagnetic phenomena. 
The problem is just to find an analytic expression for the statistical mechanical 
partition function 



where E is the total energy of the system (and the sum is over any configuration of 
the system) , from which all the thermodynamic quantities of interest (including the 
mean magnetization) can be derived. This was effectively done by Ising for magnets 
arranged into a linear chain, but he "succeeded in showing that the assumption 
of directed, sufficiently great forces between two neighboring atoms of a chain is 
not sufficient to explain ferromagnetism" [7], since no non-zero magnetization is 
predicted at any temperature. This occurs since the ordered state is unstable 
against random thermal fluctuations tending to destroy spin alignment, and Ising 
showed (incorrectly) that such negative results applies also in three dimensions. 



M [6]. 




(2) 



config. 
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This problem was not properly solved by Heisenberg in 1928, although his differ- 
ent (non constant) form of the interaction energy depends not only on the arrange- 
ment of the elementary magnets, but also on the speed with which they exchange 
their places. It was correctly addressed only in 1936 by R. Peierls ^llj who, against 
quite a commonly expressed opinion "that the solution of the three-dimensional 
problem could be reduced to that of the linear model and would led to similar re- 
sults" , finally showed that the above-mentioned statistical model of ferromagnetism 
gives non-zero spontaneous magnetization at finite temperature in two (or more) 
dimensions. 

Returning to the Heisenberg theory, his derivation of the Weiss formula ^ 
proceeded through standard calculations on the partition function 12]). The first 
step was to find the distribution of the energy levels of 2n localized electrons (one 
valence electron per atom) to first order in the coupling by evaluating the exchange 
integral between each electron and its nearest neighbors. Given the exponential 
decrease of the exchange integrals, the couplings among more distant electrons 
are, in fact, negligible. Formal group-theoretic methods were used to evaluate the 
energy levels and, in order to calculate the partition function, their fluctuations 
were neglected. This last assumption effectively meant the replacement of the 
energy levels of the spin multiplets of the microscopic magnets forming a given total 
spin S by the average energy corresponding to S. Subsequently, with a "somewhat 
arbitrary" assumption, Heisenberg introduced a Gaussian distribution of energies 
for S around its mean energy Eq with variance AEq; as he himself wrote, "it is 
of course only to be expected that the temporary theory sketched here offers but 
a qualitative scheme into which ferromagnetic phenomena will perhaps be later 
incorporated" [7]. 

Subsequent improvements, both on the mathematical formalism and on the com- 
prehension of several physical effects, came out between the end of 1920s and the 
beginning of 1930s, mainly due to F. Bloch, J. Slater and H.A. Bethe. Firstly, 
Bloch and Slater concentrated (independently) on the understanding of the physi- 
cal basis of the Heisenberg model. While considering a free electron gas (instead of 
the bound electrons contemplated by Heisenberg in his Heitler-London-like [12] ap- 
proach) for studying the role of conduction electrons in ferromagnets (for obtaining 
its thermal properties), Bloch avoided Heisenberg's Gaussian assumption and found 
a contribution coming out "when two electrons are close together, in the neighbor- 
hood of the same atom" [13] . Such an important effect (the exchange integral can 
become negative, this decreasing the possibility to have ferromagnetism), in fact, 
cannot be taken into account in the Heitler-London-like approach used by Heisen- 
berg. Furthermore Bloch replaced the difficult group-theoretic method employed 
by Heisenberg with the simpler determinantal method introduced shortly before 
by Slater in his theory of complex spectra [13] . The great advantage of writing a 
many-electron wavcfunction in the form of a Slater determinant was just to take 
into account from the very beginning the correct antisymmetry properties imposed 
by the Pauli principle, thus making sure that no two electrons in the system can 
be found in identical quantum states. This novel approach directly led Bloch to 
discover the so-called "spin waves" [TS], i.e. the states of the system corresponding 
to single or few spin flips in the fully aligned ground state, in the region of low 
temperatures, where both the Weiss phenomenological theory and the Heisenberg 
calculations couldn't be applied. 
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All these results by Heisenberg, Bloch and Slater were reviewed (and criticized) 
by W. Pauli in his comprehensive review of magnetism delivered at the October 
1930 Solvay conference in Paris [TCj. 

This meeting was also attended by E. Fermi who, very likely, was responsible 
of the spreading of the novel results about ferromagnetism within his group in 
Rome. Although Fermi and his close collaborators never worked and published on 
this argument, several people visiting the Institute of Physics in Rome in the early 
1930s did, including, primarily, H.A. Bethe [T7]. In Rome, he calculated eigenvalues 
and eigcnfunctions (at first approximation) for any value of the resulting magnetic 
moment of the whole crystal (in the onc-dimensional case of a linear chain), by 
introducing the famous Bethe ansatz for the case of two interacting spin waves. 
Other minor (to some extent) contributions came out from D.R. Inglis [I8j-[2Q], who 
visited Rome in 1932, and from a member of the extended Fermi group, G. Gentile, 
who provided [3T] a thorough justification of Bloch's results on elementary groups 
(Uebergangsgebiet) of aligned spins inside the ferromagnet [22j . Inglis, instead, 
considered the problems about the applicability of Heisenberg's approximation on 
the Gaussian distribution of energy eigenvalues [18] and, later, the effect of the 
spin-orbit interactions in ferromagnetic anisotropy [19' and the problem of the lack 
of orthogonality for the electronic wavefunctions of two neighboring atoms , 20j . 

Here we are not interested to the quite irrelevant question about the possible 
origin of such works within the Fermi group in Romejj but rather report on an 
unknown but very interesting contribution on ferromagnetism by another distin- 
guished member of that group, E. Majorana, who much contributed to the achieve- 
ments of several results in different fields of research, results that were often related 
to other people (for a recent review, see |23| and references therein). As quite usual 
for Majorana, many of these achievements were never published by him (and, in 
most cases, even his friends and colleagues were not aware of them), so that they 
are practically unknown. Fortunately enough, however, Majorana's calculations 
have been preserved in his personal notebooks [24] [25] and we are, thus, able to 
reconstruct, in particular, Majorana's contribution on ferromagnetism |25| . This is 
presented in detail in the following section, while wc comment on its relevance and 
some intriguing aspects in Section 3. Finally, our conclusions and main results are 
summarized in Section 4. 

2. Majorana theory of ferromagnetism 

The Majorana contribution begins with focusing on a set of n atoms, considered 
as magnetic dipoles (or, briefly, spins), arranged into a given geometric array at 
locations gi, (72, ■ ■ ■ 9n- Among these n elements, he assumes that a number i of them 
has a spin parallel to a given direction (spin-up), while the remaining n — i elements 
are antiparallel to that direction (spin- down). Spin-up elements are characterized 
by the quantum numbers: 

ri,r2,...ri t T • ■ • T, 
while those with spin-down are characterized by: 

ri+i,ri+2, . . .r„ ii ... i . 



Just to quote a different example, we point out that Bloch and Gentile worked together in 
Leipzig when, in 1931, they were both at the Institute headed by Heisenberg 
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The wavefunctions A(ri . . .ri\ri-^i . . .Tn) describing such a system have the form 
of a Slater determinant; in the following expression, the 5-functions (or, more ap- 
propriately, the Kronecker symbols) 5{sk T 1) guarantee that spin-up/spin-down 
elements have s = ±1 spin (in some unit), respectively: 

V'ri(gi)<5(si - 1) ... ^/'ri(gTi)<5(Sn - 1) 



A{ri . ..ri\ri+i ...Tn) 



Since there are 



V'n(ai)^(si - 1) 
'^i'r«+i(9i)<^(si + 1) 

V'r„(9l)'5(si +1) 



'4}ri+Aln)5{Sn + 1) 



(3) 



i\{n — i)\ 

ways to have i spin-up and n — i spin-down on a total of n elements (the order of 
Ti . . .ri or rj+i . . . r„ is not important), the wavefunctions describing the system 
are the following: 

■^i{l'l^1'2 ■ ■ ■ ''i ki+l ! • • • 

(4) 

^rirl, rl\rlj^-^,rl). 

Denoting with H the hamiltonian operator describing the interaction acting on each 
element, the electrostatic interaction potential Vq is given by: 



Hipi{qi)ilJi{qi)'4)2{q2)'>p2{Q2) ■ ■ ■ilJn{qn)ip„{qn) dqi . ..dqn. 



(5) 



The exchange energy between states characterized by r and s quantum numbers 
(or, according to Majorana's language, the exchange energy between r and s orbits) 
is instead: 



ki - 92! 

with Vrs = Vsr- The energy eigenvalues are, then 



{qi )i^s (<?! ) V-",. (92 )^s (92) dgidg2 



Hrr 



(6) 



(7) 



r—r\ s=r . 



(the second term describing the action of given elements on themselves, while the 
third one that of given elements on other elements), while for m ^ n: 

-Vrs, for a transition from to An by exchanging 
the opposite intrinsic orientation in the orbits 
Hrrin = i V'r and tps, (8) 

0, for the other cases. 

Majorana now assumes that (e is a constant) 

£, for neighbor atoms, 



Vrs = 



(9) 



0, 



for distant atoms, 
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each atom having a number a of neighbor atom. The energy change in the ferro- 
magnetic case is: 



E = H-Vo + Y,Vrs (10) 



2 



whose eigenvalues are: 



'1 'i+l 

e 



and, for m ^ n: 



-Vrs, 

Emn = { (12) 

0. 

Let y = y{N) be the number of configurations in which N spins can be arranged 
(at a given temperature) among each other, subject to the normahzation condition: 

N ^ 

In order to evaluate the partition function of the system, the precise determina- 
tion of the energy levels (llOp . corresponding to a given spin state of the system, 
is required, but this is virtually impossible due to the large number of elements 
involved. In his theory, Heisenberg therefore made the approximation that the dis- 
crete succession of such energy levels can be replaced by a continuum, and that such 
a continuum is distributed according to the Gaussian error law around the average 
energy for the given spin state of the system. In his notebooks, Majorana makes a 
similar assumption: "Can we consider E as diagonal, in a statistical sense? Let us 
assume that it can be." The distribution function, now subject to the normalization 
condition 

y diV = ( " j , (13) 
is then assumed to be picked around a pronounced maximum 

yo = y{No) (14) 

as in the case of a Gaussian function. This means that, with a maximum probability, 
we have yo quantities A in ([5]) corresponding to A^o- 

Then, Majorana continues: "In each of the quantities A we exchange randomly" 
(due to temperature effects) "an orbit t with an orbit J,; the quantities A change 
into quantities B: 

Ai Bi, 

(15) 

Ar > Br. 

Statistically, the set of quantities B coincides with that of the quantities A. If we 
perform the transformation (jlSp , the quantities B corresponding to the yo quantities 
A will be distributed between iVo — 2a and A^o + 2a." 
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No 2N,j 



Figure 2. A sketch of the contmuous distribution function y{N). 



Let p2r be the probabiUties that one among the mentioned B quantities corre- 
sponds to A'^ = iVo + 2r, with —a < r < a. The "average increment" in the number 
of ahgned spin is 

a 

AiVo = 51 P2r (16) 
r— — a 

and can be evaluated as follows. On average, one of the i spin-up elements has 
No/i spin-down and a — No/i spin- up elements as neighbors, while one of the n — i 
spin-down elements has Nq/ {n — i) spin-up and a — Nq/ [n — i) spin-down elements 
as neighbors. By performing the mentioned transformation (|15|) . we have 

A]Vo = 2a-2('^ + ^V (17) 

Since we have a total of £ = 2a spins around Nq, the probability that a given 
element be a spin-up is assumed by Majorana to be 



" 2 4a 
while for the probability to have a spin-down: 



' ■ (18) 



1 AiVo , , 

Thus, the probability to find m = a + r spin-up elements that alternate every 
£ — m = a — r spin-down elements is given by 



9"(l-#-", (20) 

since q is evidently a random binomial variable. By replacing in (|20p the corre- 
sponding quantities in (fTS]) . (|19p . we then get: 

1 AAoV^Vl AAoV"'' (2a)' 



^^'■-U^^J U" W {a-r)\{a + r)\ ^^^^ 
(such a law for the probability is called "normal" by Majorana). 

According to statistical mechanics, the ratio of the probabilities (jlSp and (fT9| is 
just given by the Boltzmann factor: 

^-i— =e-ra'. (22) 
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However, Majorana proceeds in a slightly different way. "Assuming that, for a 
restricted range, 

2/(iVo + 1) = 2/0 e^, 

(23) 

y{No±a)=yQ e±^^ 

the condition that y{N) does not change while we pass from quantities A to quan- 
tities B can be expressed as: 

a 

^ P2.e-2^'^ = 1, (24) 

r— — a 

that is solved by: 

2a 

which is just Eq. ([22|) . e^ being the Boltzmann factor. In the continuum limit, K 
may be replaced by y' /y. "It follows that, by considering y as a continuous function 
of TV: 

„ N (\ 1 
2 - + 



a \i n — I 



and setting 



we have: 



Q=-f- + ^)- -^ ^ (27) 



a \ I n — I 



Such a differential equation is solved by means of separation of variables, with the 
change of variable 

'=^-'- (^^) 
After few calculations, Majorana finds: 



The constant c corresponds to: 



2;(0) = y ( - ) = c (32) 



a 
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(a limit procedure is, obviously, adopted), and may be calculated from the normal- 
ization condition ()13|) . Here, in order to evaluate the integral, Majorana approxi- 
mates the function y{N) to a Gaussian curve. 



(note that y{N) takes its maximum at N = l/a, and y{l/a) — 2^/"c), obtaining 
{N' = N -1/a): 

From this, we get: 

2 

n \ I l \ " [a 



^J\V V ^ ^''^ 

The solution found is, then, worked out with a specific numerical example, for 
n = 10, i = 3, a = 4 (corresponding to a = 5/42, 2/a = 16.8 and c = 0.0002046), 
and a numerical tabulation for it is reported in the notebook. Majorana also pushes 
further the approximation (|33|) by expanding, in a Taylor series, logy rather than 
directly y, finally obtaining: 

6 15 28 45 ^ ' 

with k a given constant. 

At a given temperature T, the partition function of the system is evaluated from 
the statistical term 

y 

which now replaces the distribution function y. The Boltzmann factor may be 
written as 

e"-^ = e"& = e"-^^ (36) 

with 

The maximum of the probability distribution function (corresponding to = Nq) 
may be evaluated by taking 

^log(^ye--j^O 

(and checking, as done by Majorana, the sign of the second derivative). From (as 
above, N' = N -1/a): 

log(ye-^) =k-l-LN'--logil-a^N'')+N'log^-^, (38) 
\ y a a 1 + aN' 



we obtain: 



aNo 



l=e^ 



No = -TT71-TS (39) 
a(e^ + 1) 
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(note that for L = one re-obtains A'o — 1/ck). By setting N — Nq in ([3T|) we have 

yo = cf^Ve-^o (40) 



and 

yoe-S=c f^Ve^^oe-^A^o^,, — ) . (41) 



More in general, by expanding Eq. ([38]) in a Taylor series up to second order terms 
in {N - No) and using ((351) (gT]) we get: 

ye-^=c * e-^"(---°)^ (42) 

From the expression (l34|) we finally obtain: 



y e ''T dN = 



n 



t \ 2 



1 + e--^\ 2%/e^ 



(43) 



(in the last passage, Eq. ([27)1 has been used). 

At this point, Majorana considers the case where a magnetic field H is applied 
to the system, and puts 

"'^ <«' 

{n denoting the magnetic moment of a given atom). The additional magnetic 
energy is ruled by the net number of spins aligned along the direction of H (n^ ) or 
counter-aligned along the same direction (n^): 

- (n| - nx)^H = -(2i - n)^H. (45) 

Note that the local magnetization, in units such that the absolute saturation value 
is one, is given by 

5*=^ ^ = , 46 

+ ni n 

so that the mean magnetization of the system, i.e. "the ratio S between the mag- 
netic moment under the influence of the field H and the saturation magnetic mo- 
ment" , is: 



i 



12 



S. ESPOSITO 



The (logarithm of the) partition function 

Z=Y. ^'(S), (48) 

where p(S) is the probabiUty to have the magnetization S, in the continuum Hniit 
may then be evaluated from 

log /e-^ e^'^(2»~")2/diV - ^ Hn~i) ^ + "'"^ +2Mi-i\ogi 
J n 2 

— {n — i) log(n — i) + constant terms, (49) 
where the Stirling approximation for the factorials (in the binomial coefficients) 

log n! ~ n log n — n 

has been used. Here the omitted terms do not depend on the number i of spin-up 
elements. The condition for having a maximum in the magnetization (in the case 
of ferromagnetism) may then be obtained "by taking the derivative of (j49p with 
respect to i and equating the result to 0" : 

2(n — 2z) 1 — e^^ 

a— log h 2M - log i + log(n - i) = 0, 

n 2 

log : = 2M + log , 

n — I n 2 

or, by using (pS)) . 

= 2^ + 2a^log — (50) 
For small H and large T, this reduces to: 

25 = 2^ + 2a51og— (51) 
kl 1 + e fci" 

"For T lower than the Curie point, for a given value of H there exist 2 values of S 

which, for not extremely high H, are practically equal and opposite." 

Let now s > be the local magnetization of the system for H = 0: 

1 + s 2 

log ^ 2as log _^ . (52) 

i — s 1 -f e 

The influence of the magnetic field H on the magnetization S may be evaluated by 
putting S ~ s + As in Eq. (|50p : at first order in As we have: 

/xH 

As = — ^ 2 . (53) 

2 " " 

1 — 1 — e fci' 



3. Some implications and further contributions 

The main body of the Majorana contribution about ferromagnetism, present in his 
Quaderni, is all that reported in the previous section (apart from some omissions 
of scratch calculations and few other interesting points; see below). Here we point 
out some peculiar features of the Majorana theory, directly obtainable from the 
calculations reported but not explicitly mentioned in the notebooks. Nevertheless, 
even though in a form probably different from that given below, it is likely that 
Majorana was aware of them. 
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Figure 3. Graphical solution of the implicit equation ([54| . ob- 
tained from the intersection of the two curves y ~ log and 
y = rjx for different values of the parameter 77. A non-vanishing 
spontaneous magnetization (corresponding to the point P) exists 
only for rj > 2, that is a > 2 and T < T^. 



Let us start with stressing that, besides the trivial solution 5 = 0, Eq. (|50p ad- 
mits also another non-vanishing solution for H = 0, depending on the temperature. 
In fact, the equation: 



log ^ =VS, 



rj — 2a log ■ 



1-5 ' ' ■ "1-He-w 

has a 5 ^ solution if the following condition is fulfilled: 

2 



2a log ■ 



— > 2, 



e < 2e 



1. 



1 -I- e kT 

This condition is non trivially satisfied only if 2e^» — 1 > 0, that is 

1.44. 



1 



a > 



log 2 



In such a case, the condition ([55)1 gives: 



T < 



log (2c- i 



(54) 



(55) 



(56) 



(57) 



Thus, if the temperature is lower than the critical value given in Eq. ([57)) . the 
model studied by Majorana predicts a spontaneous (H = 0) magnetization 5 7^ 0, 
i.e. ferromagnetism, exactly as in the phenomenological Weiss theory, although in 
a form different from that coming out from Eq. ([T]). In addition to this, and as 
a necessary condition for it to occur, the inequality must be fulfilled, that is 
the integer number of nearest neighbor atoms has to be at least 2, a prediction 
not contained in the Heisenberg model (see the Introduction). Note also that, in 
the Majorana view, the quantity a effectively may count the number of physical 
dimensions (see below) of the crystal structure of the ferromagnetic body and, in 
this respect, such a result anticipates by several years what definitively proved later 
by Peierls [TT] . 
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Figure 4. Some applications, with different geometries, of the 
Majorana model of ferromagnetism. From the top: (a) geometric 
arrangement of n = 30 elementary magnets (with a = 1 near- 
est neighbor) into a linear chain; (b) triangular chain of n = 24 
elementary magnets having a — 2 nearest neighbors; (c) possi- 
ble geometric arrangements of n = 6 elementary magnets, each of 
them with a = 3 nearest neighbors. 

For some unknown reason, Majorana did not study further his model of ferro- 
magnetism, probably regarding it incomplete]^ but he nevertheless considered some 
applications to ferromagnetic materials with different geometries, corresponding to 
different numbers i of aligned spins on a total of n and different numbers a of nearest 
neighbors. Here we only mention these peculiarities, leaving aside less interesting 
numerical calculations of several quantities appearing in the model [25]. 

The simplest example is that of a linear chain but, contrarily to what usually 
considered (with a = 2), Majorana used a = 1 as depicted in Fig. H^. Explicit 
numerical calculations for this case referred to n = 60 and i = 10. 

The two dimensional case, with a = 2, was instead implemented by means of a 
triangular geometry, as showed in Fig. |4Jd, with n = 24 and i = 6. 

Finally Majorana also considered the case with a = 3 (and rt 6, i = 3), realized 
with two possible, different geometries as illustrated in Fig. 

While it is unquestionable that such examples are the simplest ones in the Ma- 
jorana model, it is unknown to what particular materials Majorana thought when 
discussing them. Note, in fact, that the only elements known at that time to exhibit 
ferromagnetic properties were iron, cobalt and nickel, whose crystal structures are 
apparently different from those considered by Majorana (body centered cubic for 
Fe, face centered cubic for Fe and Ni, hexagonal for Co). However, the realization 

^Contrary to several other topics, none about ferromagnetism was reported in the Volumetti 
I24| , this being an indication of the fact that the author did not regard his calculations as conclusive 
(see the Preface of Ref. ,24,). 
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of unusual geometric arrangements, even in usual materials, could be well explained 
in terms of different interactions responsible of the effective equilibrium positions 
of the elementary magnets. In this respect it is quite illuminating the case of a 
"triatomic molecule forming an equilateral triangle" in a state analogous to that of 
ferromagnetism, studied by Inglis [20] and, especially, that of pyrrhotite (an iron 
sulfide mineral with a variable iron content). The same author modelled this el- 
ement in terms of "one d-electron (or several independent d-electrons) per atom, 
with fairly large spin-orbit interaction (but not as large as the interatomic exchange 
energies), in a hexagonal crystal which has the iron atoms in successive principal 
planes not directly opposite one another, but staggered, or else far apart" [19]. If 
the reasoning of Inglis is really similar to that of his friend Majorana, the model 
prepared by this scholar was probably intended for providing a suitable framework 
for explaining peculiar phenomena (such as ferromagnetic anisotropy and others) 
not considered in the bare Heisenberg model. However, since no further certain 
indications exist in the Majorana notebooks, we prefer to not lucubrate on this, 
but rather stop here our analysis. 

4. Conclusion 

In this paper we have set forth and analyzed a theory of ferromagnetism developed 
by Majorana as early as the beginnings of 1930s, but never published by him. Such 
a theory differs substantially (in the methods employed and in the results obtained) 
from the well-known Heisenberg theory of 1928, although the general framework is 
the same. In fact, in both cases, the guidelines are those of the statistical Ising model 
of 1925, and the forces responsible of ferromagnetism are just derived from the 
quantum mechanical exchange interactions. Differently from Heisenberg, however, 
Majorana described the system of spin- up/spin-down elementary magnets in terms 
of wavefunctions written in the form of a Slater determinant (as did by Bloch and 
by Slater in 1929), in order to take directly into account the constraints imposed 
by the Pauli exclusion principle. 

The distribution function y{N) for the configurations of the system of spins 
modelling the ferromagnet is, then, explicitly evaluated in the continuum limit 
(that is, for a continuous distribution of the energy levels of the system, as in 
the Heisenberg model). Indeed, from the observation that given configurations of 
spin-up/spin-down elementary magnets are statistically equivalent to those where 
a spin-up is replaced by a spin-down (or viceversa) due to thermal fluctuations, 
Majorana calculates the probability for those configurations of spins and, by solving 
a differential equation, he obtains the expression for y{N), from which the partition 
function of the system may be deduced. It is also interesting to mention that 
Majorana gives a series expansion of the distribution function obtained which goes 
well beyond the usual Gaussian approximation (with terms up to order 10). 

The introduction of an external magnetic field acting on the elementary magnets 
allows the calculation of the appropriate partition function of the system and, from 
a maximum problem (on the probability for realizing a given magnetization) , a non- 
linear equation for the mean magnetization S of the system is derived. Although 
such equation is different from that deduced by Heisenberg (and that deduced early 
by Weiss) , it nevertheless admits a non trivial solution S even in the absence of 
an external magnetic field, provided that the temperature of the body is lower than 
a critical (non-vanishing) value. The main features of ferromagnetism are, thus, 
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obtained in a more transparent manner with respect to Heisenberg and subsequent 
(for instance, Bloch) derivations. 

In addition, and differently from these previous models, a condition on the num- 
ber of nearest neighbors directly follows in the Majorana theory (instead, in the 
Heisenberg theory, for example, a condition on this number is only supposed on 
the basis of available experimental data), anticipating in some respects a later key 
result by Peierls. 

Finally, some applications of the theory to peculiar arrangements of the ele- 
mentary magnets (linear chain, triangular chain, etc.) are as well considered by 
Majorana. 

From what presently discussed, then, it comes out clear the relevance of the 
contribution by Majorana to the theory of fcrromagnetism which, as happened for 
several other contributions by this scholar in different fields of physics research |23j , 
unfortunately remained unknown until now. This evidently urges to look into the 
several thousands of pages of his personal research notes [24] [25], whose careful 
study has been undertaken only recently and which will certainly deserve, once 
more, further interesting results. 
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